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The supranova model for 7 -ray bursts (GRBs) has recently gained popularity. In this scenario the 
GRB occurs weeks to years after a supernova explosion, and is located inside a pulsar wind bubble 
(PWB). High energy protons from the PWB can interact with photons from the rich radiation field 
inside the PWB or collide with cold protons from the supernova remnant, producing pions which 
decay into ~ 10 — 10® TeV neutrinos. The predicted neutrino flux from the PWBs that host the 
GRBs should be easily detectable by planned 1 km® detectors. 
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The leading models for Gamma-Ray Bursts (GRBs) 
involve a relativistic outflow emanating from a compact 
source. The 7 -rays are attributed to synchrotron (and 
possibly also inverse-Gompton) emission from relativis¬ 
tic electrons accelerated in internal shocks within the 
ejecta. The ultimate energy source is rapid accretion 
onto a newly formed stellar mass black hole. Long du¬ 
ration (> 2 s) GRBs, which include all GRBs with ob¬ 
served afterglow emission and ~ 2/3 of all GRBs, are 
widely assumed to originate from a massive star progen¬ 
itor. This is supported by growing evidence for GRBs 
occurring in star forming regions within their host galax¬ 
ies Q . The leading model for long duration GRBs is the 
collapsar model [ 3 , where a massive star promptly col¬ 
lapses into a black hole, and forms a relativistic jet that 
penetrates through the stellar envelope and produces the 
GRB. An interesting alternative, though somewhat more 
debated, model is the supranova model 0 , where a su¬ 
pernova explosion leaves behind a supra-massive neutron 
star (SMNS), of mass ~ 2.5 — 3Mq, which loses its rota¬ 
tional energy on a time scale, tgd ~ weeks to years, and 
collapses to a black hole, triggering the GRB. The most 
natural mechanism by which the SMNS can lose its rota¬ 
tional energy is through a strong pulsar type wind, which 
is expected to create a pulsar wind bubble (PWB) Hi 

In this Letter we consider the neutrino production in 
the PWB in the context of the supranova model for 
GRBs. These have ~ 10 TeV energies and are emit¬ 
ted over the time tsd between the supernova and GRB 
events, for tsd ^ 1 yrj while for tgd ^ 20 yr, that is 
required for GRBs with typical afterglows (for spherical 
PWBs) i, they are emitted mainly over the first few 
years after the supernova. For tsd ^ 0.1 yr or for GRBs 
not pointed towards us, the iz’s would not be accompa¬ 
nied by a detectable GRB. We find that the ~ 10 TeV 
neutrino fluence from a PWB at z ~ 0.1 — 1 implies 
~ 0.1 — 10 upward moving muons in a km^ detector. We 
expect ^ 10 nearby {z < 0.1) PWBs per year. The neu¬ 
trino signal from an individual PWB for which the GRB 
jet is pointed towards us, and is therefore detectable in 
7 -rays, is above the atmospheric neutrino background. 


Neutrino production in the PWB. — The main mech¬ 
anisms for neutrino production are ( 1 ) photomeson in¬ 
teractions between relativistic protons and nonthermal 
photons inside the PWB (py), (2) p-p collisions between 
the PWB protons and cold SNR protons (pp), (3) p-p 
collisions between shock accelerated SNR protons and 
cold SNR protons (s/i); p-p collisions between two rela¬ 
tivistic protons inside the PWB are unimportant. The 
pulsar wind is expected to have a comparable energy in 
protons and in pairs. The proton velocities are ran¬ 
domized at the wind termination shock, and are expected 
to be quasi-monoenergetic, with a random Lorentz fac¬ 
tor 7 p = 10 "^-® 7 p 7.5 (i.e. an energy Cp = 307 p, 4,5 TeV) 
since this is the value needed to explain GRB afterglow 
observations with tsd ^ 20 yr (for a spherical PWB) H- 
However, 7 p can be larger for tsd = 0.1tsd,-i yr < 20 yr. 

The fraction, Xp: of PWB protons that can reach the 
SNR and may interact with SNR protons, is uncertain. 
The ratio of the proton Larmor radius, Rl, and the SNR 
radius, Rb, is ^ 4 ■ ^ 1; where 

is the fraction of PWB energy in mag¬ 
netic fields, and the protons should largely follow the 
magnetic field lines. Polarization maps of the Crab neb¬ 
ula 0 suggest that the magnetic field in plerions is or¬ 
dered on scales ~ Rb. This may suggest a reasonably 
high Xp- Plasma instabilities can reduce Xp, while a ra¬ 
dial flow along the poles (jets) can increase Xp- Therefore, 
the value of Xp can range anywhere from ~ 1/2 to <C 1 . 

The average optical depth for p-p collisions between 
a proton from the PWB, that propagates in the radial 
direction, and the protons of the SNR shell is 


rad _ O-pp^SNR _ . r, f ^SNR ^ ^_2 ,-2 

PP AnRlrUp ' V 10 Mg ; "^.-i 


( 1 ) 


where Msnr and fib = 0.ll3b,-i are the mass and the 
velocity (in units of c) of the SNR shell, and cjpp « 50 mb. 
The typical optical depth to photomeson interactions is 
Tp.y « O.lA'jp l , SO that Tpj > for 7 p > 10®. The 
rotational energy lost by the SMNS prior to its collapse, 
Arot = 10 ^ 3^53 erg, is comparable to its total rotational 
energy 0. Typically AJsnr ^ 10 Mg, implying fib ^ 
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0.1 01 . Even though initially Pb is probably ^0.1, the 
large pressure inside the PWB accelerates the SNR to a 
velocity such that its kinetic energy becomes ~ E^ot 0 
0 . During the acceleratiom Rayleigh-Taylor instabilities 
are expected to develop 0 Q that may condense the 
SNR shell into clumps or filaments, as is seen in the Crab 
nebula. This increases the interface between the PWB 
and the now fragmented SNR, which helps increase Xp- 
The interaction with the pulsar wind drives a colli¬ 
sionless shock into the SNR shell 0 , that crosses it on 
a time ^ tgd, and can accelerate protons up to energies 
Ep.max ~ 10^® eV, with a roughly flat ep{dN/d€p) energy 
spectrum. These protons can collide with cold SNR pro¬ 
tons, providing an additional channel for neutrino pro¬ 
duction. The neutrino energy spectrum will follow that 
of the protons, where « ep/20. The total energy in 
these protons is a fraction /p ~ 0.1 — 0.5 of Erot- 
The pion luminosities in the different channels are 
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where = ( 2 / 3 )^p_ 2/3 is the fraction of wind energy in 
protons, and (Z^^) is the fraction of the total proton 
energy lost to pion production through p-p collisions (p -7 
interactions). A fraction rjp « 0.2 of the proton energy is 
lost to pion production in a single collision (interaction). 
The energy loss in multiple collisions can be approxi¬ 
mated by a continuous process, where ep = Cp^ exp(—r/). 
On average, n = — ln(l — r]p)~^ « 4.5 collisions produce 
one e-folding, and r/ = Tijn ~ TijA.'b, where i = pp, pj. 
This implies — exp{—T-) « 1 — exp(—0.22ri). 

The shock going into the SNR shell can produce a tan¬ 
gled magnetic field with a strength similar to that inside 
the PWB. Thus, once a PWB proton enters the SNR 
shell, it is likely to stay there over a time > fgd- This 
can increase the effective value of Tpp by up to a fac¬ 
tor of ~ [Pbi^R/Rb)]~^ ~ 100, compared to (Eq. 
P). Therefore, Tpp > 1 and ~ 1 for tsd ^ 1 yr, 
Tp-^ > 1 and - 1 for Ed < O.Ody^ f^ yr, and Tpp oc 
Tpj oc We have Tp-y ^ Tpp for Xp ^ 10® — lO’’. 

When Tpj > 1, the PWB protons lose all their energy 
via photomeson interactions before they can reach Rb 
and collide with the SNR protons. When Tp^ 1 , LPJ> ^ 
LP^, since the protons lose only ~ 20 % of their energy 
in a single py interaction, and in half of the cases are 
converted to neutrons, that can reach the SNR (since 
they are not effected by magnetic fields, and as long as 
7 p > Rb/cTn ~ Pbtsd/Tn, where « 900s is the neutron 
mean lifetime) carrying ~ 1/2 of the PWB proton energy, 
and collide with the protons there (typically Tpp > Tp^). 

So far, all the parameters were estimated at the time of 
the GRB, t = tsd- However, since Ti oc t~°‘ where app = 
3.5 {Tpp oc t/R^) and ap.y = 2.2 0 , even if Ti{tsd) < 1, at 
early times, t <ti = tsdTi(4d)^/“, Ti > 1 and // ~ 1. The 
mean value of /* over 0 < t < Gd is ~ ri(tsd)^^“a/(a — 1 ), 
and most of the neutrinos are emitted at 0 < t < 


The pions decay on a very short time scale, before they 
can suffer significant energy losses. The pion energy is 
divided roughly 1 : 1:1 ( 1 : 2 : 1 ) between 7 -rays 

and other products, for p-p (p- 7 ) reactions. The initial 
flavor ratio 4)^^ : is 1 : 2 : 0 , however, due to 

neutrino oscillations we expect a ratio of 1 : 1 : 1 at the 
Earth. The VpVp fluence over the lifetime of the PWB is 

_ QfS: _ e/.t S.o x erg 

?p/r Xpipff^ fpf^^ dl,,/{l + z) cm 2 ’ ^ ^ 

where z and d^ = 10^®<^^28 cm are the cosmological red- 
shift and luminosity distance of the PWB, respectively. 
The average flux is just the fluence divided by tsdi 3'iicl the 
fluence in 7 -rays is fPP = 2fPP, /f = 2 //^% 

The average Thomson optical depth of the SNR shell is 
(tt) = {<7t/< 7pp)Tlf- « 13t™® « (tsd/0.4yr)“2. Thereby 
fPP ~ 1 requires Tpp ~ lOOr™® > 1 and {tt) ^ 0 . 1 , while 
1 requires (tt) > 1007“4 f. If > 1 > low energy 
photons cannot escape. High energy photons, above the 
Klein-Nishina limit {t-y > rrieC^) have a reduced cross 
section and can escape if > TrmeC^ « MeV. 

Since e.y is typically < uieC^ in the prompt GRB, it would 
not be detectable for a uniform SNR with Gd ^ 0.4 yr. 
The synchrotron self-Gompton GRB emission could still 
escape. However, the optical depth of high energy GRB 
photons to pair production with the low energy PWB 
photons is T^-y > 1 for e-y > MeV 0 . This leaves 

a narrow range of photon energies that can escape the 
PWB, e-y ~ 5 — 30 MeV, for the parameter range where 
neutrino production is most efficient (Gd.-i ^ 1 )- H the 
SNR shell is clumpy, then tt ^ {tt) at the under-dense 
regions, relaxing the constraints related to tt- Moreover, 
fPP ~ 1 for tsd ^ 1 yr, so that efficient neutrino produc¬ 
tion via p-p coflisions can occur while all photons with 
£7 ^ 3(tsd/lyr)^ GeV can escape for 1 < tsd,-i ^ 10 (or 
4 < tsd ,-1 ^ 10 for a uniform shell). 

Detecting the GRB emission helps distinguish between 
i/'s from the PWB and i/'s from the atmospheric back¬ 
ground, since the PWB ^'s arrive from the same direction 
and are correlated in time with the GRB photons. 

The 7 r®’s that are produced decay into e-y^ ~ Cp/lO ~ 
37 p, 4.5 TeV photons. These high energy photons will 
produce pairs with the low energy PWB photons. The 
high energy pairs will, in turn, upscatter low en¬ 
ergy photons to high energies, etc., thus creating a pair 
cascade. This enables some photons to escape the sys¬ 
tem, even if initially Tyy{ey^o) > 1 , after having over- 
gone ^ log 2 [T 7 .y(e-y,o)] scatterings, and shifted down to 
an energy Cy ^ eyfi/Tyy{eyfi), such that Tyy{Rb) < 1 
and Tt min(l, TOeC^/e 7 ) ^ 2. Since Tyy{Rb) < 1 for ty > 
90tgd,-i MeV, then for 2 < tt < 180tgd,-i the escape of 
the photons is limited by Tyy{Rb), and the photons that 
escape have a typical energy Cy ~ MeV. These 

photons should accompany the high energy neutrinos. 
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When the SMNS collapses to a black hole and triggers 
the GRB, the pulsar wind stops abruptly. However, the 
relativistic protons in the PWB or SNR can still produce 
pions, and the neutrino emission via p-p collisions would 
not stop abruptly after the GRB, but would rather decay 
over the expansion time, ~ tsd, over which the hot pro¬ 
tons lose most of their energy via adiabatic cooling, and 
Tpp decreases. Neutrino emission via photomeson inter¬ 
actions typically decays on a shorter time scale since the 
PWB radiation field is produced by fast cooling electrons. 

The neutrino spectrum (for each of the 3 v flavors) for 
a mono-energetic proton distribution is given by 0 


$ 1 . 

K 


dN, 1 

de,ydAdt Aci, 




Ul/2)LJp' 

\m)L7, 
I 2.7/r ' 
\3.6xpfr, 


(1 + zy 


cp^L 


(4) 


g^(x)dxj (5) 


2/3^53(1 + 2)^ 10 


0^p,4.5^sd, —® 


where gTr(x) = (1 — x)^'^ + exp(—18a;)/1.34, and Cy 
ranges between m.^(? « 25 MeV (or ~ 7cm?7i7rC^/4 ^ 

GeV for p-p collisions) and ep/4 « 1 A^p^d.h TeV. 
The normalizations for py and for pp are given separately. 
The original spectrum of the 7 -rays produced in the 7 r° 
decay should be roughly similar to that of the neutrinos. 

The probability of detecting a muon neutrino in a ter¬ 
restrial ice detector is ss 1.3 x 10“®(eiy/TeV)^, with 
/3 = 2 for e,y < 1 TeV and /3 = 1 for e,y > 1 TeV M- 
For Cj/ > 10^® eV, P^p oc The expected number of 

neutrino events during the life time of the PWB is 


Np = A mm[{l + z)tsd,T] j ^uS^u)Pyii.{(^y)de^ , {&) 

NP^ = 0.05/rCp.2/3i?53(l + ^)dZ28(4l/km2) , (7) 

NP/ = 0.07/rXpCp.2/3i?53(l + ^)rfZ28(4l/km2) , (8) 

N/ = 0.05fPPfpE,/l + z)dl^/A/km^) , (9) 

where A and T ^ 1 yr are the area and integration time of 
the detector, and we have assumed T > (l-l-z)tsd- These 
expected numbers of events are valid when absorption 
in the Earth is not important, and apply to GRBs that 
point at the detector horizontally, and to < 100 TeV 
neutrinos from GRBs at the hemisphere opposite of the 
detector. Absorption may reduce the number of ev > 
100 TeV events, which are expected only for 7 ^ > 10®. 

Implications. The atmospheric neutrino background 
flux is ~ 10“^(eiy/TeV)“^ ® cm“^s“^sr“^. The an¬ 
gular resolution, 0, of the planned neutrino telescopes 
IGEGUBE and NEMO will be 0.7° [TJ and 0.3° 
respectively. The number of background events above 
a neutrino energy e,y, in the range ~ 1 — 10® TeV, 
is iVbkg(> e.) - lO-2(0/O.3°)2(e,/lOTeV)-i ®(t/O.lyr) 
per angular resolution element, over a time t. There is 



FIG. 1: The diffuse neutrino flux from different sources is 
compared to the sensitivities of present and future telescopes. 
The expected flux from PWBs that host GRBs, due to the pp 
and P 7 mechanisms (see text for details), is shown for 3 differ¬ 
ent distributions of 7 p among the different PWBs {thin solid 
lines)-. dPr/dlog( 7 p) oc 7 “, 10"' < 7 p < lO’’, a. = —0.5, 0, 0.5. 
The dotted lines show the Waxman-Bahcall cosmic-ray limit. 


better signal to noise ratio for larger Cp = 'jpUipc'^, which 
are possible for the parameter range where neutrino emis¬ 
sion is most efficient (Gd ^ 1 yr). Thus, the neutrino sig¬ 
nal from a single PWB with Gd 0.1 yr, at z ^ 1, could 
be detected with high significance above the background. 

GRBs that occur inside spherical PWBs with Gd ~ 
0.1 — 1 yr would have a peculiar and short lived after¬ 
glow emission (such events may be related to X-ray 
flashes, for which no afterglow emission was detected 
so far). Since they occur for a wide range of tsd val¬ 
ues, we expect their rate to be similar to that of typical 
GRBs (i.e. 10® yr“^ that are beamed towards us and 

^10® yr“® that are beamed away from us). Therefore, 
a km^ neutrino detector should detect ~ 100 yr"'^ neu¬ 
trinos correlated with GRBs. Furthermore, ~ 10 yr-i 
close by {z < 0.1) PWBs are expected, where each PWB 
will produce at least several events, and would thereby 
be above the atmospheric background, even if the GRB 
cannot be detected in 7 -rays, as it is beamed away from 
us or since tt > 1 for Gd ^ 0.1 yr. In the limit of very 
small Gd (;$ a few hours), this reduces to the chocked 
GRBs in the collapsar model • 

The diffuse muon neutrino flux from PWBs that host 
GRBs is shown in Fig. 1, for the P 7 -I- pp mechanisms, 
using the flux normalization from Eq. © with fr = 
fPP = 2xp = 1, and a rate of 10® yr ® at z 1. It 
is slightly below the existing u ppe r bound on the diffuse 
flux established by AMANDA Il4l| , and may be detected 
by AMANDA II. It is above the Waxman-Bahcall (WB) 
cosmic ray limit This is not a problem as these 

are obscured sources. Since the diffuse flux is the sum 
of the contributions from all PWBs, its energy spectrum 
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reflects the distribution of 7 p among the different PWBs 
(see Fig. 1), and could teach us about this distribution. 
The spectrum is typically quite different form that of 
other sources, which would make it easier to identify if 
it is detected. The diffuse flux from the sh mechanism 
is ~ 10“^/p GeV cm“^ s“^ sr“^, and extends over 

~ 10® — 10^® eV. The exact value of the diffuse flux for 
our 3 neutrino production mechanisms (py, pp and sh), 
is somewhat uncertain; the largest uncertainty is for pp. 

Our modeling of the neutrino emission from PWBs 
that host GRBs also applies to normal PWBs, like the 
Crab or Vela, when they are very young and have not 
spun down considerably, t ^ tsd- The main differences 
are that ifrot ~ 10®^ erg, Pb ~ 0.01, tgd ~ 100 yr and 
Rb oc t. Since there are ~ 100 times more regular PWBs, 
the total energy budget is comparable. Photomeson in¬ 
teractions become unimportant, but Opp = 2 and ~ 1 
at t < tpp ^ 20 yr so that time averaged over 
0 < t < tsd (f^^) 2(tpp/tsd) 0.4, and neutrino 
production via p-p collisions is rather efficient. The dif¬ 
fuse flux from normal PWBs is therefore expected to be 
comparable, to that from PWBs that host GRBs. The 
neutrino emission from normal PWBs was recently cal¬ 
culated assuming larger proton energies, resulting in 
a diffuse neutrino spectrum extending to higher energies 
but with similar e1{dN/de^) peak flux levels. Since nor¬ 
mal PWBs are much more common, galactic PWBs like 
the Crab and Vela might be detected by planned km^ 
detectors even though t > Rd for these sources M 

In order for GRBs to have typical afterglow emission, 
we need tgd ^ 20 yr, 7 p < 10® (for spherical PWBs) 0, 
implying fP'^ < fPP, (fPP) < 0.25 and Np ~ 0.01 km“^ 
events per PWB. The expected detection rate of VpDp 
from all PWBs hosting GRBs with typical afterglows, 
that are detectable in y-rays (i.e. pointed towards us), is 
~ 10 km“^ yr~^. However, most of these u’s will arrive 
~ 20(1 -I- z) yr before the GRB, at 0 < t < 4(1 -|- z) yr. 
Around the time of the GRB Jpp ~ 10“® and the rate of 
events is Np ~ 10“"^ km“^ yr“^. Since the z/’s from the 
PWB are emitted isotropically, the diffuse neutrino flux 
from these sources, ^ 2 x 10“^ GeV cm“^ s“^ sr“^, 
is dominated by GRBs that point away from us. This 
diffuse flux is somewhat lower than that expected for 
PWBs with tsd ^ 1 yr and normal PWBs, and is below 
the atmospheric neutrino background (since y^ < 10®). 

The expected number of events per PWB with tsd ^ 
1 yr at z ~ 1 is ~ 0.1 km“^, with neutrino energies 
El/ ~ 1 — 10® TeV. These z/’s are emitted over a time 
tsd ~ 0.01 — 1 yr before the GRB, and the emission de¬ 
cays over a similar time after the GRB. This emission 
is therefore easily distinguishable from emitted either 
simultaneously with the GRB or < 100 s before the 
GRB 0 , as predicted for the collapsar model of GRBs. 
The number of neutrino events from the PWB is at least 


an order of magnitude larger than from the prompt GRB, 
and is much larger than from the afterglow |2(l|. There¬ 
fore, this GRB precursor neutrino signal is one of the 
best candidates for an e arly detection with the planned 
km^ neutrino telescopes |2]J. The detection of a neutrino 
signal from a PWB can serve to distinguish between the 
different progenitor models for GRBs. 
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